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Abstract 

We propose a new model describing the production and the establishment of the stable 
gradient of the Bicoid protein along the antero-posterior axis of the embryo of Drosophila. 
In this model, we consider that bicoid mRNA diffuses along the antero-posterior axis of the 
embryo and the protein is produced in the ribosomes localized near the syncytial nuclei. 
Bicoid protein stays localized near the syncytial nuclei as observed in experiments. We 
calibrate the parameters of the mathematical model with experimental data taken during 
the cleavage stages 11 to 14 of the developing embryo of Drosophila. We obtain good 
agreement between the experimental and the model gradients, with relative errors in the 
range 5 — 8%. The inferred diffusion coefficient of bicoid mRNA is in the range 4.6 x 
10^^^ — 1.5 X 10^^^ m^s~^, in agreement with the theoretical predictions and experimental 
measurements for the diffusion of macromolecules in the cytoplasm. We show that the 
model based on the mRNA diffusion hypothesis is consistent with the known observational 
data, supporting the recent experimental findings of the gradient of bicoid mRNA in 
Drosophila [Spirov et al. (2009) Development 136:605-614]. 



Introduction 



In Drosophila early development, bicoid mRNA of maternal origin is deposited in one of 
the poles of the egg, determining the anterior tip of the embryo, (Frigerio et al. . 1986 



Driever and Niisslein- Volhard , 1988). The deposition of the mRNA is done during oo- 
genesis and is transported into the oocyte along microtubules, (Saxton, 2001). After 



fertilization and deposition of the egg, and during the first 14 nuclear divisions of the 
developing embryo, bicoid mRNA of maternal origin is translated into protein in the 
ribosomes. 

During the interphases following the 11th nuclear division up to the 14th, the concen- 
tration of Bicoid protein distributes non-uniformly along the antero-posterior axis of the 
syncytial blastoderm. Bicoid has higher concentration near the anterior pole of the em- 
bryo, and its local concentration decreases as the distance to the anterior pole increases. 



This is called the Bicoid protein gradient, ( Niisslein- Volhard 1992). 

As, during oogenesis, bicoid mRNA is deposited near the anterior pole of the em- 
bryo, it is implicitly assumed that Bicoid protein is produced in the ribosomes of the 
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nuclei localized near the anterior pole of the embryo, and then diffuses through the syn- 
cytial blastoderm. Eventually, this protein diffusion could be facilitated by the absence 
of cellular membranes in the syncytial phase of the developing embryo. Driever and 
Niisslein-Volhard (Driever and Niisslein- Volhard , 1988) argued that the protein gradient 
is generated by protein diffusion and degradation throughout the embryo. Later, Niisslein- 
Volhard (Niisslein-Volhard, 1992, 2006) emphasized that the Bicoid protein diffuses away 
from the site of its production, the local mRNA deposition region. The theoretical pos- 
sibility of this mechanism of morphogenesis goes back to the work of Alan Turing in the 
fifties (Turing, 1952), and has been further discussed and analyzed by Wolpert ( jWolpertj 



1969), Crick (Crick, 1970) and Meinhardt (Meinhardt, 1977). Experimental measurements 



within mammalian cells (Wojcieszyn et al. 



1981 Mastro et al. , 1984) and theoretical anal- 



ysis (Crick, 1970) suggested that diffusion coefficients of macromolecules in the cytoplasm 
are in the range 10~^^ 



10 



-13 j^2g 



However, there are several open questions related with the establishment of the stable 
gradient of protein Bicoid along the antero-posterior axis of the embryo of Drosophila. Ex- 
perimental observations during cleavage stages 11 — 14, and before cellularization that oc- 
curs at the end of cleavage stage 14, show that Bicoid protein is always localized around the 
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This fact suggests that dispersal effects driven by molecular collisions with Bicoid protein 
— Brownian motion — do not play a significant role in the establishment of the gradient 



of Bicoid. As Brownian motion is the driven mechanism of diffusion dispersal, (Murray 
1993, chap. 9), it is difficult to understand how the diffusion of a protein produces a 
strongly localized protein concentration around the syncytial nuclei and during successive 
mitotic cycles. On the other hand, as Bicoid protein is produced and attains a steady 
state during the first cleavage cycles, its localization near the nuclear envelops suggest 
that ribosomes are also localized near the nucleus. If ribosomes were not localized near 
the nuclear envelopes of the syncytial nucleus, protein in the inter-nuclear regions of the 
embryo would be observed. 

As argued by Kerszberg and Wolpert (Kerszberg and Wolpert, 2007), there is not 



a clear experimental evidence of protein degradation, a necessary mechanism for the 
establishment of a steady protein gradient in models based on protein diffusion. 

Houchmandzadeh et al. ( Houchmandzadeh et al.. 2005), reported the constant Bicoid 
protein concentration during cleavage cycles 12-14, suggesting the stability of protein con- 
centration during an important developmental period. On the other hand, the protein 
diffusion hypothesis lead to some quantitative contradictory facts. For example, in re- 
cent experiments, the hypothetical inferred cytoplasmic diffusion coefficient of the Bicoid 
protein during the cleavage stage 13 is of the order of 0.3 x 10~^^ m^s~^, (Gregor et al. 



2007). However, during the first cleavage stages of the developing embryo, the Bicoid pro- 



tein reaches a steady state in 90 minutes (end of cleavage stage 9), (Gregor et al., 2007), 
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and a simple estimate with the Houchmandzadeh et al. ( Houchmandzadeh et ai. 2005) 



model shows that the diffusion coefficients must be of the order of 2 x 10 



-12 j^2g 



(Gre- 



gor et al, 2007). This value for the diffusion coefficient is one order of magnitude larger 



than the value inferred from experiments. This discrepancy between model estimates and 
observation needs a clear explanation, (Reinitz, 2007). 

Here, with a mathematical model, we show that the observed gradient of the Bicoid 
protein can be explained by the diffusion of bicoid mRNA, and Bicoid protein stays local- 
ized near the nuclei of the syncytial blastoderm of the embryo of Drosophila. This explains 
the absence or the very low level of Bicoid concentration in the regions between the nu- 
clei during the first stage of development of Drosophila. We determine a scaling relation 
between the mRNA diffusion coefficient, the embryo length and the mRNA degradation 
rate, enabling the precise determination of the diffusion coefficient of bicoid mRNA. In 
this model, it is not necessary to introduce the morphogen degradation hypothesis for 
protein, and the steady gradient of protein is reached after the complete translation of 
mRNA of maternal origin. 

The mRNA localization mechanism in the embryo of Drosophila has been analyzed 
experimentally by several authors, and (Saxton, 2001) argues that the relative small size 
of mRNA suggests that random diffusion and specific anchoring to the cytoskeleton in a 
target area might suffice for localization in the syncytial blastoderm. Cha et al. (Cha et 



ai, 2001) reported rapid saltatory movements in injected bicoid mRNA in the embryo, 
followed by dispersion without localization. Other effects of diffusing mRNA has been 



reported by Forrest and Gavis (Forrest and Gavis, 2003) for the nanos mRNA. More 
recently, Spirov et al. (Spirov et al.. 2009) have shown that a bicoid mRNA gradient 



exists along the antero-posterior axis of the embryo of Drosophila., completely changing 
our current views of this Drosophila developmental pathway. The model presented here 
corroborates these experimental facts, is consistent with the experimental facts and ob- 
servations, and fits the experimental data with high accuracy. 



Results 

We now derive a mRNA diffusion model and we show that experimental protein gradi- 
ents are well fitted in this framework. This shows that a mechanism of mRNA mobility 
(diffusion) is enough to explain protein gradients. 

A mRNA diffusion model 

It is an experimental fact that mRNA of maternal origin is deposited in a small region 



of the embryo of Drosophila., defining the anterior pole of the fertilized egg, (Driever and 



Niisslein- Volhard , 1988). After the deposition of bicoid mRNA, we assume that bicoid 



mRNA of maternal origin disperses within the embryo and this process is simultaneous 
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with the successive cleavage stages. Then, the protein is produced in the ribosomes that 
are near the nuclear membranes of the nuclei in the syncytium. Representing by R{x, t) 
the concentration of bicoid mRNA along the one-dimensional antero-posterior axis (x) of 
the embryo, and by B{x,t) the concentration of Bicoid protein, the equations describing 
the production of Bicoid from mRNA are, 

= -dR + D^ 

dx^ (1) 

= aR 

where a is the rate of production of Bicoid from mRNA, d is the degradation rate of bicoid 
mRNA, and D is the diffusion coefficient of bicoid mRNA in the cytoplasm. Using the 
mass action law, this simple model is straightforwardly derived from the rate mechanisms, 

R^B + R 

and then the diffusion term is added to the mRNA rate equation. In general, one molecule 
of mRNA can produce more than one molecule of protein, implying that d < a. If one 
molecule of mRNA produces one molecule of protein then, in the mean, we have d = a. 

We consider that the length of the antero-posterior axis of the embryo is L, and so 
X E [0,L]. We take zero flux boundary conditions, §f(a; = 0,t) = §f(a; = L,t) = 0, 
and If (x = 0,)!:) = §f = L,t) = 0, for every t > 0. The protein initial condition is 
B{x,t = 0) = 0, and the initial distribution of mRNA is, 

A>0 if 0<ii<x<i2<L 
R(x,t = 0) = { (3) 
otherwise 

where A, ii and £2 are constants. The function R{x, t = 0) describes the initial distribution 
of bicoid mRNA of maternal origin deposited in the region of the embryo [ii, £2] C [0, L]. 
The concentration of mRNA of maternal origin deposited in the embryo is then A{i2 — £i)- 
In this model, bicoid mRNA has a fixed initial concentration, and the Bicoid protein does 
not degrade. 

Equation ([T]) with the initial condition ([s]), and the zero fiux boundary conditions 
define the mRNA diffusion model. This model is linear, and has solutions that can be 
determined explicitly. Now, we will show that, within this simple model, Bicoid protein 
attains a gradient like steady state along the embryo. 



By standard Fourier analysis techniques, see for example (Alves and Dilao, 2006) or 
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(Haberman, 1998), the solution of the first equation in ([T| is, 



R{x,t) 



A- 



+2-4 E 



COS 



i=l 



nn 



sm 



mT£2 



sm 



rnrii 



(4) 



The solution of the second equation in ([T]) is, 

B{x,t) = B{x,t = 0) + a f R{x,s)ds 

Jo 



(5) 



In the limit t —>■ oo, the equilibrium or steady solution of the Bicoid protein is calculated 
from Q and ([s]), and we obtain. 



Beq{x) 



L 



cos 



1=1 ai 



sm 



sm 



(6) 



where, 



A- , 02 = Ct- 



and we have introduced into ([s]) the protein initial condition B{x,t = 0) = 0. A simple 
calculation shows that the solution (|6| can be written as, (Alves and Dilao, 2006), 



Beq{x) 



o2a2/L 



ai X, 
- cosh(a2-, 



smh(a2— j — smh(a2 — 



where. 



I(x) 



ai (e-'^2(£i-a;)/L 



-a2{£2—x)/L 



-a2(x—ii)/L 



ai (e-'^2(x-£2)/L 



)/2, if x<ii 

if £i < a; < 
if X > £2 



-a2{t.2—x)/L 



(9) 



g-a2(a;-^i)/Lj ^2 

The steady solution (|8])-([9]) describes the gradient of the Bicoid protein along the 
antero-posterior axis of the embryo of Drosophila. This solution depends on the set of five 
parameters ai, 02, L, ii and £2, to be calibrated with experimental data. The constants 
ai and given by ([7| define scaling relations of the embryo. 

In order to compare the model predictions with the experimental data, the next step 
is to calibrate the parameters of the mRNA diffusion model ([Sj-Q with the available 
experimental data for the gradient of the Bicoid protein. 
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Calibration of the mRNA diffusion model with the experimental 
data 

To calibrate the parameters of the mRNA diffusion m odel ([8|)-([9|) with the experimental 



data, we use the data available in the FlyEx database, (Poustelnikova et al. . 2004 Pisarev 



et al. . 2009). We considered Bicoid gradients for cleavage stages 11-14, and we make the 



additional assumption that during these cleavage stages, the concentration of the Bicoid 



protein is in the steady state, ( Houchmandzadeh et al.. 2005). 



We have fitted the data sets of the FlyEx database with the equilibrium distribution 
of Bicoid protein given by The fitted funct ions are represented in Figure pi 

In this figure, we show the concentration of Bicoid protein along the antero-posterior 
axis of Drosophila, for several embryos and in consecutive developmental stages. Due 
to the particular form of model prediction (|8|-([9]), the fitted parameter values are a = 
{ai,a2,ii/L,i2/L)- The parameter values of the different data sets are shown in Table [l| 
and correspond to the global minima of the fitness function x^l^^), introduced below in 



(15) 



Table 1: Fitted model parameters for the protein Bicoid antero-posterior dis- 
tributions 





ai 


^2 


h/L 


h/L 


a/y2 /52 
V Am/ ^max 


n 


d (s-i) 




Aa{i2 - h)/L 


a) 


abl8 (11) 


345.2 


4.69 


0.03 


0.20 


0.06 


30 


8.8 X 10- 


-4 


5.2 X 10-2 


b) 


abl7 (12) 


894.4 


4.50 


0.06 


0.14 


0.06 


70 


8.1 X 10' 


-4 


5.8 X 10-2 


c) 


abl6 (13) 


684.2 


5.51 


0.06 


0.15 


0.08 


152 


1.2 X 10" 


-3 


7.4 X 10-2 


d) 


abl2 (14-1) 


927.6 


4.82 


0.06 


0.14 


0.08 


309 


9.2 X 10" 


-4 


6.8 X 10-2 


e) 


abl4 (14-2) 


3414.7 


4.38 


0.05 


0.07 


0.08 


314 


7.7 X 10- 


-4 


5.3 X 10-2 


f) 


ab9 (14-3) 


1191.6 


4.39 


0.05 


0.10 


0.07 


343 


7.7 X 10- 


-4 


4.6 X 10-2 


g) 


adl3 (14-4) 


470.4 


3.02 


0.06 


0.19 


0.05 


324 


3.6 X 10- 


-4 


2.2 X 10-2 


h) 


ab8 (14-5) 


3271.7 


4.25 


0.08 


0.09 


0.07 


332 


7.2 X 10- 


-4 


2.4 X 10-2 



Parameters values that best fit the experimental distribution of Bicoid protein shown in 
Figure |2] with the equilibrium distribution ([8])-(|9|. In the first column, we show the data 
sets and the corresponding cleavage stages. The parameters ai, 02, ii/L and £2/L have 
been determined with the swarm algorithm described in the Materials and Methods 
section. The lengths of the embryos have been rescaled to the value L = 1. a/ Xm/ ^max 
is an estimate of the relative error of the fits, and n is the number of data points in the 
corresponding graphs in Figure |2] The parameter d has been determined by ([T]) with the 

of bicoid mRNA and L = 0.5 x 10-^ m. 



10 



-11 j^2g 



estimated diffusion coefficient D 
The parameter ^40(^2 ~ £i)/L has been determined with (10). 



The quality of the fits of Figure [2j has been evaluated with the fitness function (15). 
Denoting by B^ax the maximum value of each experimental data set, the mean relative 
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error of a fit is estimated by the quantity, \/Xm/ ^maxj where Xm = niin^g^ In 
Table [T| we show the mean relative errors of the fits, and the number of points (n) in 
each data set. The mean relative errors between the theoretical predictions (|8])-([9]) and 
the experimental data sets of Figure [2] are in the range 5% — 8%, showing a remarkable 
agreement between the model prediction and the experimental data. 

To determine the values of the parameter d in Tab le [1} we have fixed the embryo 



length to the value L = 0.5 x 10~^ m, ( Niisslein- Volhard , 2006, cap. iv). For the diffusion 
coefficient of bicoid mRNA, we have chosen the value, D = 10~^^ m^s~^, as estimated 
below in (14). By Q, (i = alD/L'^, and the value of the degradation rate d depends on 
the choices made for D and L. 

As the experimental data is given in arbitrary light intensity units, the initial value of 
the bicoid mRNA concentration is also arbitrary. However, it is plausible to assume that 
the total amount of initial bicoid mRNA deposited in the embryo does not change too 
much for different embryos. So, in order to estimate the total amount of bicoid mRNA in 
the embryos, using the first relation in ([t]), we have calculated the quantity, 

Aa{e2 - ii)/L = aid{i2 - ii)/L, (10) 

where A{i2 —ii) is the total amount of initial bicoid mRNA deposited in the embryo, and 
the rate a should not change too much for different embryos. Therefore, if the quantity in 



(10), does not change too much for different data sets, it is an indication of the ability of 



the model to describe data sets with different phenotypes. In fact, as shown in Table [T| 



the quantity (10) is almost constant among embryos, even if {£2 — ii)/L shows a large 
variability, as is the case of the fits in Figure [2j 

From the fitted values shown in Table [T| each of the calculated parameter values 
Aa{i2 — C-ij/L, 0,2 and d have the same order of magnitude for the different cleavage 
stages. In fact, the parameters defined in the kinetic mechanisms ([2]) are independent 
of the cleavage stage and, therefore, their values must depend only of the phenotypic 
characteristics of the analyzed embryo. The similarities between the parameters Aa{l2 — 
£i)/L, 02 and d for different cleavage stages is an indication of the consistency of the 
theoretical model proposed here. 



Determination of the diffusion coefficient of bicoid mRNA 

To determine the value of the diffusion coefficient of the bicoid mRNA, we use the in- 
formation that Bicoid protein reaches a steady state in approximately T seconds. So, 
we integrate the two equations in ([T]) along the embryo length, and using the zero flux 
boundary conditions, we obtain. 
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where R and B are the total amount of mRNA and protein in the embryo, respectively. 
The differential equations (jTTl) have the solutions. 




R{0)e 

a 
d 



-dt 



R{0) (1 - e"''*) 



(12) 



Assuming that the steady state of the Bicoid protein is attained after T seconds of devel- 
opment, and that 95% of the mRNA has been translated into protein, by (12), we have 
the development time relation, B(T)/{aR{0)/d) = 0.95 = (l — e~'^^). From the previous 
relation we obtain, d = — log(0.05)/T. Therefore, by ([T]), the diffusion coefficient is. 



D = d 



log(0.05) L2 



(13) 



As Bicoid protein attains the steady state at the end of cleavage stage 9, in approximately 
T ~ 90 X 60 seconds, (Gregor et al. , 2007), with the data in Table [l| we have 02 G [3, 5.5], 
and with the choice L = 0.5 x 10~^ m, by (13), the diffusion coefficient of bicoid mRNA 
is in the range. 



D e [4.6 X 10"'^ 1.5 X 10""] (14) 

These estimates, as well as the numerical fits of Figure |2| are consistent with the 
theoretical predictions for the order of magnitude of the diffusion coefficients of large 



molecules {bicoid mRNA) in the cytoplasm (Wojcieszyn et al., 1981 Mastro et al, 1984). 



Discussion 

We have proposed a new model describing the production and the establishment of the 
stable gradient of the Bicoid protein along the antero-posterior axis of the embryo of 
Drosophila. In this model, bicoid mRNA diffuses along the antero-posterior axis of the 
embryo and Bicoid protein is produced and stays localized near the syncytial nuclei as 
observed in experiments. 

We have calculated the steady state of the Bicoid protein along the antero-posterior 
axis of the embryo of Drosophila, and we have calibrated the parameters of the mRNA 
diffusion model with experimental data taken during cleavage stages 11-14. After the 
calibration of the model with experimental data, we have predicted the initial localization 
in the embryo of the bicoid mRNA of maternal origin (parameters ii and £2 in Table 1), the 
Bicoid protein concentration profiles along the embryo, and the bicoid mRNA degradation 
rates. The mean relative errors between the theoretical prediction of the Bicoid protein 
steady state and the experimental data (Figure [2]) are in the range 5% — 8%, suggesting 
the effective validation of the model proposed here. 

A simple estimate gives a diffusion coefficients D for bicoid mRNA in the interval 
[4.6 X 10~^^, 1.5 X 10~^^] m^s~^. This estimate is calculated with the parameters found 
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in the calibration of the experimental data, with the model prediction formulas for the 
steady state, and with the additional assumption that the gradient of Bicoid protein 
is reached at the end of cleavage stage 9, (Gregor et al., 2007). The determination of 



the diffusion coefficient is strongly dependent of the time duration of the cleavage cycles 
and therefore, it has a large error that is difficult to quantify. On the other hand, as 
the steady state solution of this bicoid mRNA diffusion model depends on the scaling 
parameter 02 = \/ dL'^/D and 02 is determined by fitting the steady states of the protein 
profiles, the experimental determination of the degradation rate of bicoid mRNA leads to 
a more precise estimate of the diffusion coefficient. 

The calibration and validation of the mRNA diffusion model shows that the mech- 
anism of establishment of the gradient of Bicoid protein observed in Drosophila early 
development can be justified by a diffusion hypothesis for mRNAs. The mathematical 
model considers that bicoid mRNA diffuses along the embryo and the translated protein 
stays localized near the syncytial nuclei, as observed in Figure [2j In this model, pro- 
tein degradation is not considered and proteins do not diffuse along the embryo. The 
model proposed here explains the experimental data for protein gradients, and the com- 
mon assumption that morphogen gradients are obtained with a balanced and continuous 
production and degradation of proteins is not necessary. The low level of Bicoid concen- 
tration in the intranuclear regions of the embryo is easily explained through the ribosome 
localization near the syncytial nuclei. 

Random motion of bicoid mRNA has been observed, (Cha et al., 2001 Saxton, 2001), 



and the bicoid mRNA gradient has been recently found by Spirov et al. (Spirov et al. 



2009). The mechanism of mRNA diffusion proposed here is uni-dimensional and together 



with the very good agreement between the experimental data and the model predictions, 
we can raise the hypothesis of the existence of a mechanism of constrained mRNA diffusion 
along a network of nonpolar microtubules. This hypothesis has been discussed by Spirov 
et al. ( Spirov et al. , 2009 ) and is consistent with a mechanism based on mRNA diffusion 
along microtubules. This justifies the very good agreement between one-dimensional 
diffusion models and the observed experimental data (the motion of mRNA is observed 
along the embryo wall). Two and three-dimensional diffusion models would predict protein 
concentrations in the interior of the embryo, which is not observed. 



Materials and Methods 

To fit the sets of data points of Figure [2} we consider that each data set is approximated by 
a function B{x; a), with x G [0, 1] and where a = (cti, . . . , am) is the set of m parameters 
to be determined. We assume that the parameter space S = {a : 00 < rrii < at < Mi < 
00, i = 1, . . . m} is a compact subset of R"^. For each fixed value of the vector parameter 
a, and experimental data points {(xj, -Bexp(a^j))}r=i) consider the fitness function (sum 
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of the mean squared deviations), 



n 

i=l 



1 " 

-J2iBixf,d)-B,,,,{xi)Y (15) 



The set of parameter values that best fits the experimental data is determined from the 
global minimization condition, 



minx^fo?) 



In order to search the global minimum of the function x^('3), with a G S*, we take a 
set of p vectors dk, with k = 1, . . . ,p, randomly equidistributed in the set S. Then, for 



each dk, we compute the fitness function (15) 



To search for the global minimum of the mean squared deviation x i'^)^ '^6 swarm 
the set of p vectors dk in the parameter space S. For each vector d = (ai, . . . ,am), we 
construct a new vector d' according to the swarm rule, 

= «i + At{Mi - nii) (2^-1) (16) 

where At is a time parameter, (Mj — rrii) is a scaling constant, and ^ is a random variable 
uniformly distributed in the interval [0,1]. Then, we recalculate the new value of the 
fitness function, x^(tt')- If X^i^i') < X^('5), the parameter value d is updated to the new 
value d'. If x^(a') > X^('5); update is done. We repeat this procedure for all the 
parameter values in the search space S. 

After iterating the swarm algorithm M times for all the population of parameter 
values, we order the parameter vectors according to their fitness values, and we discard 
half of the parameters that have the worst fitnesses. We repeat this procedure s times. 
The parameter values that best fit the experimental data are the ones that corresponds 
to the minimum of X^(a). 

This simple algorithm relies on the assumption that the initial number of random 
points are equidistributed in S and they form a sufficiently dense set in S. 

In the cases in Figure [2| the convergence of the swarm algorithm for the determination 



of the global minimum of the fitness function (15) as a function of the parameters has 
been checked by graphical methods in two-dimensional sections of the parameter space 
S. In all the fits in Figure [2| the best convergence has been obtained with the swarm 
parameters, p = 1024, M = 500, s = 5 and At = 0.01, and good convergence for the 
global minima has been obtained. 

A further extension of this calibration technique using evolutionary algorithms for the 



Bicoid-Caudal protein regulation in Drosophila has been developed in (Dilao et al. , 2009). 
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Figure 1: Distribution of Bicoid protein (in blue) in the embryo of Drosophila, 
in the interphase following the cleavage stages 11 (a) and 12 (b). The images are 



from the FlyEx datasets abl8 (a) and abl7 (b), (jKozlov et al., 


2000 Myasnikova et al. 


1999 


2001 


Poustelnikova et al., 2004 


Pisarev et al. 


2009). Note the absence of Bicoid 



protein in the inter-nuclear regions of the cytoplasm. The localization of Bicoid protein 
near the nuclear envelopes suggest that ribosomes are also localized near the nucleus. 
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Figure 2: Concentration of Bicoid protein (BCD) along the antero-posterior 
axis (x) of Drosophila, in the interphase following the consecutive cleavage 
stages (cl st) 11, 12, 13 and 14. Data (connected dots) in the figures are from the 
FlyEx datasets: a) abl8 (11); b) abl7 (12); c) abl6 (13); d) abl2 (14A-1); e) abl4 
(14A-2); f) ab9 (14A-3); g) adl3 (14 A-4); h) ab8 (14A-5) and the numbers inside the 
parenthesis refer to the cle avage stage, ([Kozlov etoT 2000; Myasnikova et ai, 1999, 2001 



Poustelnikova et ai, 2004 Pisarev et a/. 2009). In the upper right corner of the figures, 



we show the corresponding gradients of Bicoid protein in the two-dimensional projections 
of the embryo. The data points of the Bicoid gradients are taken from a region centered 
around the central antero-posterior axis of the embryo. The transversal length of this 
region is equal to 10% of the maximal length of the dorso-ventral direction. The lengths 
of the embryos have been rescaled to the value L = 1. The thick black lines are the best 
fits of the experimental data with the theoretical prediction ([8])-(|9]). The gray regions 
show the initial localization of bicoid mRNA, and are defined by the fitted values ii/L 
and £2/1^- The parameters of the fits are shown in Table [l} 



